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SUMMARY 


Modeling the human as a controller of slowly responding systems 
with preview is considered. Along with control tasks, discrete non-control 
tasks occur at irregular intervals. In multitask situations such as these, 
it has been observed that humans tend to apply piecewise constant controls. 
It is believed that the magnitude of controls and the durations for which 
they remain constant are dependent directly on the system bandwidth* preview 
distance, complexity of the trajectory to be followed, and nature of the 
non-control tasks. A simple heuristic model of human control behavior in 
this situation is presented. The results of a simulation study, whose 
purpose was determination of the sensitivity of the model to its parameters, 
are discussed. 



INTRODUCTION 

Although successful operation of an airliner is now possible from 
take-off to touchdown with minimum involvement of the human pilot [1] he 
must still perform various routine checks in the course of a normal flight. 
In addition, even when flying on autopilot, constant monitoring of various 
instruments is necessary to detect any out of tolerance signals and abnormal 
occurences of any events. Further, malfunctions or changes in atmospheric 
conditions, for example, might require that the pilot take over control and 
nake course changes that are different from the preplanned trajectory. Thus, 
despite advances in automation, human control of aircraft is certainly still 
of interest. 

When the human is controlling a plant* it has been observed that 
the controls applied are not always continuous. Continuous controls are 
necessary and are observed when the time constants involved are rather small 
and the deviations from some reference trajectory must be kept within some 
close tolerance. But when the time constants are relatively large, it is 
unnecessary and also difficult to apply the right amount of continuous 
control. For slowly responding processes it is often sufficient and 
desirable to apply step-like controls intermittently. This gives an 
opportunity to observe the actual behavior of the system, compare it with the 

* This research was supported by the National Aeronautics and Space 
Administration under NASA-Ames Grant NSG— 2 119* 

75 





predicted behavior, and take corrective action. This usually prevails in a 
tracking situation where a certain length of the future command trajectory is 
available, along with the present required position. Further, applying the 
step-like controls also frees the human to engage in non-control tasks. In 
fact, this kind of Dehavior is common in process control situations and also 
has been observed in simulations of a flight management situation [2], [3]. 

When preview of the command trajectory for a certain distance 
into the future is available, it is likely that the human would apply 
step-like controls so as to minimize the future trajectory deviations rather 
than instantaneous deviations. A model which appears reasonable is one which 
updates the expected deviations of the cost over the length of the previewed 
trajectory and uses this .information along with the knowledge that it "costs” 
to change control values. The cost to change control reflects the fact that 
non -control tasks must be attended to, though they may not be of primary 
importance. The "cost” is thus due to the feeling that the non-control tasks 
would “suffer" if attention is focused away from them and on the primary task 
alone. This cost may manifest itself as a tolerance threshold for error 
below which no action is taken. A measure for the cost of not attending to 
the subsystem tasks is available as a function of various probabilities and 
costs for delay of subsystem tasks [M]. 


BACKGROUND 

Of the available models for manual control, the optimal control 
model would appear to be a suitable candidate. However, this model assumes 
that control remains non-zero at all times whereas in an intermittent control 
situation, control is zero for a significant portion of the time. Hence the 
mean fraction of time devoted solely to control, corresponding to non-zero 
control intervals, cannot be calculated with the optimal control model. For 
a given fraction of attention, the conventional optimal control model 
predicts only RMS errors and RMS control actions. While recent versions of 
the model do yield a measure of attention that should optimally be used for 
monitoring subsystems that dynamically relate to the control task, subsystem 
t*.sks that only remotely relate to the aircraft's dynamic response cannot be 
considered. Further, in multitask situations the optimal control model's 
performance criterion, which minimizes mean squared deviations, may not be 
appropriate. Finally, these approaches do not yield any predictions of the 
split of attention between control and non-control tasks or about the 
probability that the human is involved in the control of a continuous system 
at any particular instant. 

The human in multitask situations has been modeled by Walden and 
Rouse [3] as a ‘server' in a queue where 'customers' are the control and 
non -control tasks. The customers are assumed to arrive for service with 
exponentially distributed inter-arrival times (Poisson arrivals. ) Service 
times are Hrlang-k distributed. Some customers have a higher priority over 
others (e.g., control tasks over non-control tasks.) There are a total of N 
customers in the population (total number of possible tasks the human may be 
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called upon to attend), and N epacee are available J In »• «“?"« 
worst all the N systems may require service simultaneously.) hi 
can be modeled as a (M/E k /1 :PBP/N/N) queue. ( s ee references [5] or t6]for 
details about the notation. ) The queueing model predicts the Action of time 
soent in each type of task (i.e., server utilization). The emphasis in thi 
model is on the subsystem task performance. The control task is modeled in 
the sense that performing it consumes time. However, measures charac enzmg 
control performance (i.e., RMS errors) are not available. 


Ah INITIAL MODEL 


Some success has been achieved using a heuristic model to 
describe control of an aircraft (with simplified dynamics) in a horizontal 
plane. Initial computer simulations indicate that this could be a fr it 
approach. A plece-wis. atraight line nap vaa created using unifornly 
distributed random variables for the length of straight line segm n • 
magnitudes for angle of turn between segments were chosen from nine values 
( 10*-90 # ) with equal probability. The direction was chosen randomly. This 
type of Ip was designed because of the flexibility in determining the 
parameters. It is a simple matter to change the probability dist ^ b ^ 1 °"g i ? 
various parameters of the path, so that different conditions could be easily 
tested. It was assumed that the aircraft would be ^ing forward with 
constant speed. A point moved along the map corresponding to ^he tesired 
aircraft position. A distance equivalent to two time constants ahead of the 
desired position on the map was shown as preview. Only lateral motion was 
considered. Control in the horizontal plane was achieved through use of th 
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The deviation from the desired position was constantly monitored. 
A perfect internal model of the aircraft was assumed. At every time instant, 
the model calculated the future positions for the entire preview length 
assuming that the roll angle would be zero. At any point, the error was 
calculated as the deviation from the commanded position, instead of the 
perpendicular distance to the map 

It appears reasonable to expect that the importance given by the 
human operator to deviations from the desired trajectory will vary along the 
length of the trajectory. Deviations at the current time cannot be corrected 
or changed significantly due to the slow response of the system. Also, the 
expected deviations near the end of the previewed command trajectory need not 
be considered immediately, since enough time will be available in the future 
to correct these. Further, any changes "beyond the horizon" that would come 
into view soon can reasonably be ignored. Accordingly, the human might 
weight the mid-portion of the previewed command more than either end. So 
the weighting function for errors would increase to a maximum (from zero), 
about one time constant frcm the current position, and decrease again to a 
near zero value at the end of the trajectory. 



TRAOttTOMES . REFERENCE HOOEL 

Fig. I Simulation Rraults : Heuristic Mode l 


Weighted error is squared and summed over the preview leneth. If 
this predicted error function exceeds a certain threshold, ’aileron '"is held 
at a maximum value until maximum bank angle is reached. If the error is 
within the threshold, the bank angle is made zero. Constant weights on 
errors were uped for the simulation. The results are shown in Figures 2 and 
3 and appear reasonable. The time intervals corresponding to non-zero 




aileron action give a direct measure of fraction of time required for control 
which is proposed as a correlate of workload. Although the human must 
continuously monitor for cumulative error, workload due to this is assumed 
negligible compared to the workload involved in control where he must watch 
the effect of his actions more carefully. 



Fig. 3 Control and Roll Angle History for a Given Course Change 


Simulation experiments were conducted usii.,; i fractional 
factorial design to find out the sensitive/relevant parameters. A resolution 
VI design was used so that no main effect or two factor interaction is 
confounded with any other main effect, two factor interaction, or three 
factor interactions. The following parameters were assumed to affect 
performance: 


1. Dynamics ( t of the process), 

2 Average arrival rate of turns, 

3 Standard deviation 
of course changes 

4. Amount of preview 

5. Weighting function, 

6 Ihreshold on cumulative 
weighted error. 
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fMS errors and the fraction of time spent on control were used as performance 
measures. A constant function (rectangular) and a triangular function were 
used for weighting on errors* For other parameters, two extreme values were 
chosen, to obtain a total of 32 different conditions. Exponentially 
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distributed segment lengths and normally distributed angles 
for the map. For each condition, two replications were run. 
snown in Tables 1 and II. 


Table I 


Analysis for RMS Error 


Factor 

Effect 

Sum of 

DOF 

F' Ratio 




Squares 




1 

4-433 

314.409 

1 

131.037 

P < 0.001 

2 

-0.671 

7-206 

1 

3-003 


3 

1.477 

34.887 

1 

14.540 

P < 0.001 

4 

4-087 

267-247 

1 

111.381 

P < 0-001 

5 

-1.397 

31.219 

1 

13.011 

P < 0.001 

6 

0.622 

6.188 

1 

2-579 


12 

-0.268 

1.145 

1 

0.477 


13 

1.108 

19.653 

1 

8-191 

P < 0.01 

14 

4.174 

278.726 

1 

116.165 

P < 0.001 

15 

-0.886 

12.557 

1 

5.233 

P < 0.05 

16 

-0.250 

1 .004 

1 

0.418 


23 

-0.288 

1.325 

1 

0.552 


24 

0.049 

0.038 

1 

0.016 


25 

0.335 

1.794 

1 

0.748 


26 

1.315 

27.676 

1 

11.535 

P < 0.005 

34 

0.383 

2.350 

1 

0.979 


35 

1.295 

26.822 

1 

11.179 

P < 0.005 

36 

0.197 

0-619 

1 

0.258 


45 

-0.914 

13.364 

1 

5.570 

P < 0.025 

46 

0.258 

1.062 

1 

0.443 


56 

0.124 

0.246 

1 

0.102 


Average 

4 -192 





Error 


100.774 

42 



Total 


1150.313 

63 




1 1-Period, 2-Segment Length. 3-Angle of Turn, 
4-Preview Length, 5-Weight, 6-Threshold) 


of turn were used 
The results are 




Table II 


Analysis for Fraction of Attention 


Factor 

Effect 

Sum of 

DOF 

F Ratio 






Squares 






1 

0-231 

0.856 

1 

33-221 

P 

< 

0.001 

2 

-0-057 

0-051 

1 

1-999 




3 

0.017 

0.005 

1 

0.181 




4 

0.144 

0.331 

1 

12.837 

P 

< 

0.001 

5 

-0. 1 17 

0.218 

1 

8.468 

P 

< 

0.01 

6 

-0.139 

0-311 

1 

12.084 

P 

< 

0.005 

12 

0-004 

0.000 

1 

0.011 




13 

-0.018 

0.005 

1 

0.206 




14 

0.269 

1.158 

1 

44.963 

P 

< 

0.001 

15 

0.000 

0-000 

1 

0.000 




16 

0.004 

0.000 

1 

0.008 




23 

0.067 

0.072 

1 

2.783 




24 

0.030 

0.014 

1 

0-541 




25 

0-005 

0-000 

1 

0.017 




26 

0.060 

0.058 

1 

2-234 




34 

-0.074 

0-087 

1 

3-359 




35 

0.047 

0.035 

1 

1.374 




36 

-0.043 

0.030 

1 

1.155 




45 

0-005 

0 000 

1 

0-014 




46 

0.048 

0.037 

1 

1-435 




56 

0.090 

0. 129 

1 

5.011 

P 

< 

0-05 

Average 

0-256 







Error 


1.082 

42 





Total 


4-479 

63 






( 1-Period » 2-Segment Length. 3-Angle of Turn, 
4-Preview Length 5-weight, 6-Threshold) 


It can be seen that period, preview length and their interaction 
have the largest effect on performance. Different weighting functions also | 

affejt performance* In addition, RMS error is affected by the magnitude of I 

turns and various interactions. Fraction of time spent on control is a 

affected by the threshold. Higher threshold values reduce this fraction. a 

Though the interaction of mean segment length and threshold a 

affects the hMS error, segment length alone does not affect either of the a 

performance measures. This could be due to the constant forward speed in all 8 

cases, whereas the mean segment length is scaled by the time constant. For ? 

slower process, for a given threshold any error that may result takes a | 

longer tirce to reduce to zero. Since the 'vehicle' would stay away from the f 
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trajectory for a longer time* higher RHS error results. 

for the “S hl « her MS emr ls observed 

for a different control aotlM °°" 8 "J° vi8 "' ,nd " 1 * ht ««n 

also quite frequently rSiitTi In in °.=t . f 1 ^ the error <*•«*•• 
result- « M * resulting m increased control action. This attain 

a uit cyoL Z:»7ZT IT th « freehold. Ihue, behavior «.n,Tt„ 
subjects control sIok processes^^MV" 817 ^ hSS be * n obServed when naive 
one threshold above which control is actSatid posaibly be avoided by having 
which control is made zero. * 1 * actuated » and a lower threshold below 


CONCLUSIONS 

experimental situation* 51 for* use wittf human subim t inVOl * e d?velo P ment of an 
be included to simulate a < sub J ects * Non-control tasks will 

be used. «£lSK£?i « tills vlT'ZT*' * Simple aritho * tio tasks may 
possibility. Complexfty and the rate at wmTI ° f reSUlts af * e a 

^Spuia^d. 000 ^ 01 U8k «**«*« <*** be 

°" tth 

approach, it may be possible ! this regard is the satisfaction 

obtain a heuristic-based solution withth^ addition in J his frameWo »*, and 
related to the problem structure Wth Vm?*?, f 8 eM condlt ions 
realistic experiment will be deveionert \ h ® Se P mod ® ls availab l e * a more 
(UAT II). Along With the ? *1 the General Aviation Trainer II 

more realist" 8 °‘ ,ntro1 tMk ' the "“-oontrol tasks win be lad. 

presented, simulation’ results 1 for'a^set °of°oondfti r T. 00 " tro1 t,5k “*» 
trajectories were given. The controller na rt a ondi ti°ns describing various 
Only the threshold must be d£er£L3 El internal nod ^- 

period and the preview length were found to h*® 1 *! internittent control. The 
affecting performance. This will form th* > ^ "I?®* iroportant Parameters 
kith humans. The model will be refined to '..vT 3 »'° r pr0p03ed experiments 
these experiment, end then klllbaled ^ 8 * 000u ” t tl,e "•»«» of 

non -control tasks, to modal the overall multitask situation^' 1 " 8 "° d81 f ° P 
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